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Abstract 

The problem considered in this paper is to find when the non-central Wishart distribution, 
defined on the cone Vd of semi positive definite matrices of order d and with a real valued 
shape parameter, exists. We reduce this problem to the problem of existence of the measures 
m{n,k,d) defined on Vd and with Laplace transform (det s)""''^ exp tr(s~^w) where n is an 
integer and where w = diag(0, . . . , 0, 1, . . . , 1) has order d and rank k. We compute m(d— 1, d, d) 
and we show that neither m{d — 2, d, d) nor m(d — 2, d — 1, d) exist. This proves a conjecture 
of E. Mayerhofer. 

Keywords: non-central Wishart, zonal polynomials, Euclidean Jordan algebras. AMS clas- 
sification: primary 60H10; secondary 60B11. 

Abreviated title: Non-central Wishart distributions. 

1 Introduction 

The non-central Wishart distribution is traditionally defined as the distribution of the random 
symmetric real matrix X = YiY{ + ■ ■ ■ + YnY* wlrere Yi e M'*, i = 1, . . . ,n are independent 
Gaussian column vectors with the same non-singular covariance matrix S and respective means 
rrii, i = 1, . . . ,n not necessarily equal (here * means transposition). For s in the open cone Vd of 
positive definite symmetric matrices of order d and w — mim\ -I- • • • 4- ninm^ in the closed cone of 
semi positive definite matrices Vd one can readily derive the Laplace transform 

ip/ -tr(sJf)N _ I -tr(2s(/d+2Ss)-iu;) N 

^ '~ det{Id + 2i:s)^/^ ■ ^ ' 

It is important to note that in this formula the rank of w is < n. 

Exactly like the statistician who extends the familar chi square distribution with n degrees of 
freedom to the gamma distribution with a continuous shape parameter, one is tempted to extend 
the values that the power of det{Id + 2T,s) can take in the above formula. The question is then: 
given T, £ Vd and w £ Vd, for which values of p > does there exist a probability distribution on 
Vd for X such that for all s € T'd we have 

_tr(sJf)N _ I -tr(2s(/d+2Ss)-i™)7 (2) 

^ ' det(/d + 2Ss)P ■ ^' 

Call this hypothetic distribution for X satisfying Q a non-central Wishart distribution with pa- 
rameters (2p, w, S), or NCW{2p^ w, E) for short. This question was adressed in Letac and Massam 
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(2008) where we claimed (in Proposition 2.3) that such a distribution exists if and only if p is in 
Arf with 

,12 d - 2, 1 , 

Ad = {2,2,---,^}U[^,(X3) (3) 

without any restriction on Y, Cz Vd ot^ w (z Vd- This statement was actually not quite proven: 
we considered as obvious that if p is in the part {i, |, . . . , ^} of then NCW{2p, w, S) does 
exist. The gap in our proof was kindly pointed out to us, in a private communication, by E. 
Mayerhofer who later showed in Mayerhofer (2010) that the statement was not only unproven, but 
false. More specifically Mayerhofer shows that, if NCW{2p, w, S) exists, if d > 3 and if 2p is in 
{1, 2, . . . ,d — 2}, then rankw < 2p + 1. He even conjectures that rank w < 2p must hold. The 
aim of the present paper is to correct our mistake by giving a satisfactory necessary and sufficient 
condition of existence of NCW{2p, w, E): in Section 2, we show that if 2p is in {1, 2, . . . , d— 2} then 
rankw < 2p, thus proving that the Mayerhofer conjecture is true. These questions are delicate 
and Mayerhofer (2011) uses a stochastic process valued in the set of symmetric matrices in order 
to prove his main statement (reformulated in Proposition 2.4 below). The methods of the present 
paper are simpler and necessitate only a careful study of the support of the measures in Vd when 
these measures are obtained by convolution. The basic tool (Lemma 5.2) is the following. Let Mi, 
be the set of positive measures concentrated on the matrices of rank bofVd- It is not generally true 
that if a + < d, /i e Ma and v e Mt, then /i * j/ e Ma+b , but it is true if either fi or v is invariant 
by X I—)- uxu^^ for any orthogonal matrix u of order d. This result is the subject of Section 5. It 
is, however, not sufficient to prove Propositions 2.4 and 2.5: we need further information about 
the measure m{d — l,d,d) on Vd defined by its Laplace transform (det s)"'-''"^-'/^ exp( tr (s""'^)). 
We need to show that it has an absolutely continuous part. For this reason. Propositions 3.1 
and 4.4 give a description of m(l, 2, 2) and m{d ~ 1, d, d). We give there more details than strictly 
necessary but to describe the singular and absolutely continuous parts of m((i— 1, d, d), in Section 4, 
we make an interesting use of zonal polynomials. We are guided by an elementary study in Section 
3 of the case d — 2 which uses the Faa di Bruno formula only. We should also mention that the 
proper mathematical framework for this paper is that of Euclidean Jordan algebras rather than the 
linear spaces of real symmetric matrices and we make frequent references to Faraut and Koranyi 
(1994) (henceforth abbreviated FK). But working in that framework might have obscured our 
statements without adding any insight: the extension of our results to Euclidean Jordan algebras 
is straightforward. 

We would like to thank here E. Mayerhofer for pointing out our mistake and conjecturing the 
correct result, and J. Faraut for helping us with Lemma 4.1 below. 



2 Reduction of the problem: the measures m{2p, k, d) 

Let k be an integer such that < fc < d. We consider the diagonal matrix /(fc, d) with its first 
d—k diagonal terms equal to and the last k equal to 1. For p G we define the positive measure 
m{2p, k, d) on Vd such that for all s G Vd have 

J_e~'''^'^^m{2p,k,d){dx) = .^^j-i-^e*^'""'^^^^^". (4) 

Note that m{2p, k, d) may or may not exist. Formula ((22|) below gives m(l, 1, 1). For 2p > d — 1 
formula ([39|) gives m{2p, d, d). If k and n are integers such that < fc < n < d, formula ([7|) gives 
m{n, k, d). Finally m{d — 1, d, d) is computed in Sections 3 and 4. The paper will show that these 
examples are the only cases of existence. The following proposition links this measure m{2p, k, d) 
with our initial existence problem. 

Proposition 2.1. Let E S Vd, w £ Vd andp <E A^. Suppose that rank w = fc. Then NCW{2p, w, E) 
exists if and only if m{2p, fc, d) exists. 
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Proof. Assume that m{2p, k, d) exists and let us show that NCW{2p, w, S) exists. The proof is 
based on the following principle. Let be a positive measure on a finite dimensional real linear 
space E such that its Laplace transform L^{s) = e~^'*'^V(da;) is finite on some convex subset 
D{n) of the dual space E* with a non empty interior. Let a be a linear automorphism of E* and 
let b G E* such that Li^i_{a{h)) < oo. Then there exists a probability P{a,b) on E with Laplace 
transform Lp(^a,b){s) = Lf^{a{s + h)) / L ^{a{h)) . This probability P{a,b) is obtained in two steps: 
first take the image vidy) of ^{dx) by the map x i— > a*{x) = y where a* is the adjoint of o. Its 
Laplace transform is L^{s) — Lp(a(s)). The second step constructs P{a,b){dy) as the probability 
e-^^'y)v{dy)lLMb))- 

Let us apply this principle to the case where E — E* is the Euclidean space of real symmetric 
matrices of order d with scalar product {x, y) = tr {xy) and where /i is m{2p, k, d). Here = Vd- 
We take b = (2T,)~^ and a to be the linear transformation s i— a(s) = qsq* where q is an invertible 
matrix of order d such that 

2(2E)-iu;(2S)-^ = q-^ I{k,d){q*)-^ . (5) 
We have a*{x) — q*xq. The distribution P{a,b) is the noncentral Wishart NCW{2p,w,Y,) since 

Lfj.{a{s + b)) ^ 1 -tT(2s(Ia+2T:s)-^w) //jn 

L^(a(6)) det(/rf + 2I]s)f ' ^' 

The verification of ([5]) is easily done by a calculation of trace using tr (ab) — tr (ba) and (O: 

tr[((g*)-i(s + (2I])-i)-ig-i-(g*)-i(2S)g-i)/(fc,d)] 

= tr [{{s + (2S)-i)-i - 2J:)q-^I{k, d){q*y^] = - tr (2s(/d + 2Ss)^^w) 

The only thing left to prove is the existence of q satisfying ([5]). To see this, since the matrix 
2{2T,)-'^w{2T,)-'^ oiVd has rank k, we write 2(2T,)-^w{2J:y^ = uAu* where 

A = diag(0,...,0,A?,...,A^) 

with Ai > and where u is an orthogonal matrix of order d. Taking q — diag(l, . . . , 1, A^^, . . . , A^^) m* 
provides a solution of ([S])- 

The proof of the converse follows similar lines. □ 



Example: When < k < n < dwe can use the above principle for constructing NCW{n, 21 {k, d), Id) 
from •m{n, k, d). We take q — Id and b = Id/2. Since a is the identity we have therefore 

m{n, k, d){dx) = 2*'/2e2fegtr:./2^^^(^^ 2/(fc, d), Id){dx) (7) 



The next three propositions reformulate known facts in the langage of the measures m(2p, k, d). 

Proposition 2.2. Let n and k be integers such that < n,k < d. The measure ■m{n,k,d) exists 
for < fc < n < d. Furthermore, the measure m{d — l,d,d) exists. 

Proof. Formula ([7]) provides an explicit form of m{n,k,d). For 2p > d — 1 the probability 
NCW{2p, Id, Id) exists as proved in Letac and Massam (2008) Proposition 2.2. This implies that 

lim NCWi2p,Id,Id)^ NCW{d-l,Id,Id) 

p\{d-l)/2 

exists by considering the Laplace transforms (this crucial remark is due to Mayerhofer (2010)). 
From Proposition 2.1 we have the result. □ 

Proposition 2.3. Suppose d > 3. If m{d — 2,d — 1, d) does not exist then m(n, k, d) exists for no 
pairs (n, k) such that < n < fc < d. If m(d — 2, d, d) does not exist then m(n, d, d) exists for no 
n such that < n < d — 2. 
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Proof. Suppose that m{n,k,d) exists for some pair < n < A; < d. We define m'{dx) as the 
measure on Vd with Laplace transform 

[_e-''('^^m'{dx) = \__^t.[s~\nd-i4)-i(k4))] 

JVd (det s) 2 

Since the rank of I{d ~ 1, c?)) — /(fc, d) is equal to d — 1 — fc and less than or equal to d — n — 2 then 
m! exists by Propositions 2.1 and 2.2. Now we write the convolution 

m{n, k, d) * m' = to(c? — 2, d — 1, d) 

which contradicts the non-existence of m(d — 2, d — 1, d). Similarly, suppose that m(d — 2, d, d) does 
not exist and that there exists n such that < n < d — 2 and such that m(n, d, d) exists. Then 
m(n, d, d) * m(d — 2 — n, 0, d) = m(d — 2, d, d) also leads to a contradiction. □ 

The idea of the proof of Proposition 2.3 is essentially due to Mayerhofer (2010). Here is now his 
important main result: 

Proposition 2.4. If d > 3 the measure m(d — 2, d, d) does not exist. 
Here is our main result: 

Proposition 2.5. If d > 3 the measure m(d — 2, d — 1, d) does not exist. 

We will prove Proposition 2.5 in Section 6. In the remainder of the paper we develop the tools 
that lead us to this proof. They will also enable us to give a quick proof of Proposition 2.4. Let us 
emphasize the fact that Propositions 2.1 to 2.5 give a necessary and sufficient condition of existence 
of the distribution NCW{2p, w, E). Suppose that d > 3. Given T, ^ Vd , w £ Vd and p > then 
NCW{2p, w, S) exists if and only if the following are satisfied 

1. pe Ad defined by ©; 

2. if 2p > d — 1 then rankw is arbitrary; 

3. if 2p = n < d — 2 then rank w < n. 

Note that for d = 2 the probability NCW{2p, w, S) exists if and only if 2p > 1. 



3 Computation of m(l, 2, 2) 



In this section we compute to(1,2,2) using only calculus. We parameterize the cone V2 by the 
cone of revolution C = {{x,y, z) G M.^;x > + using the map ip from C to V2 defined by 



X + y z 
z X — y 



Note that tr b, c)(p(x, y, z)] — 2ax + 2by + 2cz. 



Proposition 3.1: If C = {{x,y,z) G M.'^;x > \/ y^ ^ z^} consider the positive measure ^ on C 
such that for (a, 6, c) £ C we have 



1 



-2ax — 2by—2cz 



^{dx, dy, dz), 



\/ a? — — JC 
that is to say such that the image of /i by is m(l, 2, 2). Then 

/i(da;, dy, dz) = s{dx, dy, dz) + f{x, y, z)lc{x, y, z)dxdydz 



(8) 
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where the singular part s is the image of the nieasure g{2y^y^ + z^)dydz on by the map {y, z) i— >■ 
{x, y, z) = {-sfW^^-, Vi z) with 

2 

g{z) — — cosh(2-v/z) 
nz 

and where for (x, y, z) €z C 

„2 „,2 ,2\fc °° 



2 ^ {x'-y'-z'f ^ 1 {2xy 

j(x,y,z) ^2^^ kl{k + iy. Z^^T{m + 2k+l) ml 



Proof: Let D — Recall the following differentiation formula. 

Faa di Bruno formula . If /(t) and g{x) are functions with enough derivatives, then 



where k — ki + ■ ■ ■ + kn and where the sum is taken on all integers kj > such that fci + 2fc2 + 
• • • + nkn = n. 

For a reference see for instance Roman (1980). We apply © to g defined by a; i— — — for 
fixed y, z and to f{t) = t". Noting that D^g = 0, we obtain 

We now recall that ior p > 1/2 we have 

(a2_52_^2)p ^ r{p)T{p-^)Jc y y ) y y ) 

The idea of the proof is to apply (jlip for p = n + | and to write 
1 ^ 2a I ^ (2a)"+i 



V- 2a "+^ 1 

> 7 — —r, X 3- (12) 

^ n + 1! (a2 _ ^2 _ c2-)n+i 



Va2 - 62 _ c2 Va^ - 62 _ c2 (n + 1)! [a^ - - c^y 

1 



\/ — b'^ — 

\n+l 



V^^An+ l)!n!r(n + |) ic y z j cixdydz 



„^p(n+l)!n!r(n+f) 7c 
Define 



A first step is to observe that for fc = 0, 1, . . . , n we have 



4(n) = ^ 2''2'-2- A_(^2 _ ^2 _ ^2)n^^^y^^_ (^3) 

Let us prove it by induction on k. It is true for fc = 0. Suppose that it is true for k < n and let us 
show that ([T^ is true for k + 1. Observe that for fixed (y, z) the root y^y^ + z^ of the polynomial 

X I— >■ (x^ — y^ — z^)^ has order n and this implies that ■§^{x'^ ~ J/^ ~ 2^)" is zero for x = y/y"^ + z"^. 
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Using this remark and integration by parts with V{x) — e and U{x) — ^fc(a;^ ~ z'^)"': we 
compute the fohowing integral: 



fjk+1 

-2ax " /2 „,2 „2\n. 



jix' -y'~zTdx (14) 



With we are in position to prove (fT5|) . We have 



Qk 



h+i{n) = 2a I e-^^^-^'y-'^^^ix'-y'-z^rclxdydz 



^-2by-2cz 



^-2by-2cz 



2ae~''^^ix^-y^-z'rdx 
^--rp^ dx" 

e'^'^^S^ix'-y'-z'rdx 



Qxk+i 



dydz 
dydz 



fik+l 



which proves (|13p . We will need (fT3|) only for k = n. The second step is to express /„+i(n) as the 
Laplace transform of a positive measure. We compute /n(n) using again an integration by parts. 
The new fact for fc = n is that the integrated part will not disappear and will provide a term for 
the singular measure s given in the statement of the theorem. This calculation of the integrated 
part will use dTU]). Taking V{x) — —e~'^°-^ and U{x) — -§-^{x'^ — — z^)"' , we have 



/„+i(n) = 2a/„(n) 



-2by-2cz 



R2 



gn 



^-2a.-2b,-2c.£_(^2 _ ^2 _ ,2^n^^^y^^ 

+nl f e-'''y-'^^[-e-''^^2xrr/^dydz 
where ([T5)) comes from ([TU)) by keeping only the term ^2 = 0. We finally obtain 



dydz 



(15) 
(16) 



/„+i(n) = / e-^''''-^''y-^^'^^-—{x^-y^-zy'dxdydz+nl / e-2«Vy'+^'-2ba-2c^(2^y2 + z'^Ydydz. 
Jc oa:"+i yj[2 

Note that ^ is 



p— ? (n + l)!n!r(n+|) 



-/„+i(ri). 



The last step of the proof is to represent the function on C \ dC defined by (a, 6, c) i~-> \i 
as a Laplace transform. Using the Gaussian integral in (|17p we obtain 



1 



Va^ - 62 - c2 



g-2a(«^+t,")-2b(u^-t,^)-4cu«^j^^^ 



2 



(17) 
(18) 



To derive (fTS]) observe that the map on {(u,w);m > 0} defined by y = ^2 _ ^ _ ^ 
bijection with 'R^; the same is true with {{u,v)\u < 0}. Furthermore dydz — 4(m^ + v'^)dudv = 



Ay'y^+z^dudv and therefore dwdu 



dydz 

4y/y^+z 



All of this leads to (UHl). 
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Now we use r(n + i) = i^," i and we consider the function 



I 1 



2 2 ^"^^ 1 2 

^ t" = — cosh(2\/t) 



(?^+ l)!r(?i+ §) TTt 



We then define the measure s{dx, dy, dz) concentrated on the boundary dC of the cone to be the 
image of the measure 



g{2^/y^^ + z^)dydz = ^ cosh {2^/^{y^ + zY^^)dydz 

7r-\/ y^ + z'^ 



(19) 



by the map {y, z) i— > (x, y, z) = + 2^, y, z). This s will be the singular part of the measure. 

We now concentrate on the absolutely continuous part. We will need the following formula, 
similar to PH)) and also obtained by the Faa di Bruno formula 



9a;" ^ ^ ^ ' (fc2 - 1)! 



fc2 = l 



X — X ) — ' _ (20) 



{n - 2k2) 



The absolutely continuous part of m(l, 2, 2) is given by the part of ([T2|) defined on C and its 
density is 



f{x,y, 



2 1 9"+^ 2_ 2_ 2xn 

rn + l)!7i!r(n+|)9a;»+i^^ ^ 



V2] 



J_ ^ 1 1 {x^ -y^- z^ f^-^ (2x)"-2''^ 



fc2! 

1 (2a;) 



(n- 2fc2)! 



y/^f-' fc!(fc + l)! ^r(m + 2fc+|) m 



□ 



Remark: The image by Lp of the measure s(dx, dy, dz) is concentrated on the set 5*1 C V2 of 
matrices of rank one. Any element of S\ can be written as u 







u* where u is an orthogonal 



matrix of 0(2) and Ai > 0. We can compute the image of s{dx,dy,dz) by the map 



X + y z 
z X — y 



\/y^~+^ + y 



z yj/^ + z^ — 



y J 



Ai = 2v/?y2 + z2. 



If At — {{x,y, z);2y/y^ + z^ < i}, then using polar coordinates y — AiCOSq;,z = Aisina with 



Jacobian equal to ^ , we have 

s{At) — / s{dx,dy,dz) 



At 



, 9{2Vy'^ + z^)dydz J / 5(Ai)AidAi. 

2^y^ + z^<t ^ Jo 



Since g{Xi) = cosh2-\/Ai , the image of s is 



cosh(2v Ai)l(o,oo)(Ai)dAi 



(21) 
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Now an important observation is the following: consider the measure m(l, 1, l){dX) on (0,cxi) 
whose Laplace transform is 

1 1/, ^ 1 r A"-^ 1 r 1 



=e 



y -r = / e-'^^ y — ^7 rrd)^ = ^ / e-«^ — cosh(2VA)dA 



As a consequence 



i(l, 1, l)(dA) = cosh(2VA) 1(0,00) (A)dA (22) 

V A 



and one observes that m(l, 1, 1) is quite close of (|21l) . To summarize this remark, the singular part of 
to(1, 2, 2) can be seen as the image of \pKX\m{\^ 1, V){d\\)®du by the map (w, Ai) w ^ 





from (0, oo) x 0(2) where du is the uniform probability on 0(2). This is the key to the generalization 
of m(l, 2, 2) to ra{d — 1, d, d) for c? > 2 done in Proposition 4.4 below. 

4 Computation of the measure m{d — l,d,d) 

Before stating Proposition 4.4 which describes m{d— l,d,d) we have to fix some notations, recall a 
few facts about zonal functions and polynomials and prove three lemmas. The Lebesgue measure 
dx on the space of symmetric matrices of order d has the normalization associated to the Euclidean 
structure given by (x, y) — tr (xy). Let £d be the set of sequences k — (mi, . . . , m^) of d integers 
such that nil > 7712 > ■ . . > rrid > 0. li k, € £d we consider the two zonal polynomials ^k{x) = 
^mi,...,md{x) and Ck,{x) = CK.{Id)^K{x) where C^ild) is defined below in ([M)) . To introduce 
we consider 

A,{x) = Ai(x)™i-"^A2(x)™^-™^ ...Ad_i(x)™'^-i-™''Ad(x)'"'^ 

where for x = {xij)i<ij<d a- real symmetric matrix, A.k{x) = det(a;y )i<i.j<fc. The function is 
defined by 

$K(a;) = / A^{uxu*)du (23) 



where du is the Haar probability on the orthogonal group 0{d). When x & Vd definition (1231) 
makes sense even when mi, . . . ,md are complex numbers. In that case ^mi,...,maix) is no longer 
a polynomial and is called a zonal function. To give the value of the constant C^ild) we need the 
notations £{k) — max{j; m.j > 0}, |k| — mi + . . . + m^ and 

Td{zi, ...,Zd)^W_ V{zj —) 

defined for Zj — > 0, j = 1, . . . , d. If p is a real number, we use the convention 

Vd{z + p) = rrf(2:i + p, . . . , + p)- 

If K e and if p > (d— l)/2 we define the Pochhammer symbol (p)^ = Y diK. + v) 1^ dij>) ■ If G f d 
the constant CK,{ld) is given in Muirhead page 237 formula (38) by 

ILli (2n^i +«(k) - «)' 

We never consider Ck.{x) if k ^ Ed- The exact value of CK,{Id) will be important in the proof of 
Proposition 4.4 when we shall need the formula (3) page 259 of Muirhead (1983): 



E ^^«(-) (25) 
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We are indebted to Jacques Faraut for the next lemma: 



Lemma 4.1: If x = 



Xi Xi2 
X21 X2 

numbers mi, ... , we have 



S Vd we denote [x]i = x\ £ Vd-i- Then for ah complex 

$mi,...,mrf(a;) = (detcc)™" / 'i>,n,,...,m^_d[uXU*]l)du. 



Proof: Consider v = 



vi 
1 



e 0{d) where vi e 0{d - 1). Observe that 
[vyv*]i = vi[y]ivl 



We write 



(deta;)™'' / A^,^,„^^^_^{[uxu*]i)du 



= (deta;)™'' / A„n,...,md-i{[vuxu*v*]i)du 



(26) 

(27) 
(28) 
(29) 



(det x) I '—^nii .nid 

/0(d) 



(dctx)"'' 



0(d) \ JO(rf-l) 



-^{vi[uxu*v*]ivl)du (30) 
A.m^^....,md-i{vi[uxu*v*]ivl)dvi^ du (31) 

(32) 



= (dctx)™'' / <^mu-,md-A[uXU*]i)du 
JO(d) 

In this list ([?7)) comes from the definition of ^i^{x), ([^5]) separates the roles of [uxu*]i and 
det(ii2:u*) = detx in the definition of Ai^{uxu*), ((29|) uses the fact that du is the Haar mea- 
sure, (150)1 applies (|26p to y = uxu* , (PT| uses the fact that the Haar measure dui of 0)(d — 1) has 
mass 1, (15^ comes from the definition of $mi....,md_i (2^)- Q 

Lemma 4.2: If a; G T'd then ^,ni,...,md{x^^) = $_,„^,... -m^ (x) 

Proof: Define p G 0(d) by 



P = 








1 

1 



1 







and define A* 



^d-i,rnA^) = ^rm,....,7nd-i,n.Apxp*)- We Can Write 



A. 



-umAuX ^U*)du 



A* 



A, 



^ {u*xu)du 



■md,---, — fn 



^ {uxu^)du 



(33) 
(34) 



In this list (p3|) comes from the FK formula Proposition VII. 1.5 (ii) page 127 which says 

Ami,...,md^i,md{x ) — '^-md,...,-mi 



9 



and (j34l) conies from the invariance of the Haar measure du on 0{d) by m i— )■ pu* .O 
Lemma 4.3: li x e Vd then $,„i^..._„^(a;)(det x)p = ^^^+p^,,,^^^+p{x). 
Proof: We have from the definition ^i^{x){dct x)p = 

/ Ajni,...,ma{uXU*){detuXU*)Pdu^ / A„^^+p^,,,^,na+p{uXU*)du^ ^rni+p,...,ma+p{x) O 
JO(d) Jo(d) 



Proposition 4.4. Define the singular measure r{dt) on Vd and concentrated on Sd-i as the image 
of the product measure 

(Trdetx)^/^ , , , , , 

— to(c;- l,(i- l)(dx) i^du 

1 (d/zj 

by the map from Vd^i x 0((i) to Vd which is {x,u) t-^ t — u 

C^{t) 



X 




u = uxu 



Define 



/.W = (detr^(^ - 



where £' d ^ {n ^ £d ] "m-d > 0}. Then 

m{d- l,d,d){dt) = r{dt) + fd{t)lj,Jt)dt. 

Proof. The fonction fd{t) is a well defined analytic function around x — since from the definition 
([23| of $«; the polynomial Cmi_...,m<i(a;) is divisible by (detx)™''. Therefore (det a;)~^CK(a;) is a 
polynomial when k € f^. Recall the basic fact (see FK, Lemma XL2.3 or Theorem XIV. 3.1): 

e- (det x)P'^ ^"^^^ , dx = $,(s"i)(det s)-p. (35) 

, Td{K+p) 



Note that the choice of the suitable Lebesgue measure dx is crucial in ((35|) . From ([251) we know 
that, for 2p > d — 1, the Laplace transform of m{2p, d, d) is 

/ e-*^(^")m(2p,d,d)(dx) = (dets)-?' V (36) 



Observe that the Laplace transform of fd{t)\-p^{t)dt is easily deduced from p5p and is 



e 
•Pd 



--(-)/.(^)dt (det.)-^ 



In ([55)) take 2p = d — 1. From the Laplace transform (I37p our aim is therefore to prove that the 
Laplace transform of r{dt) is 



e 
•Pd 



--(-V(do = (det.)-^i =(det.)-- E ^^tPl 



,K;e£d\£'d / Vftefd 



To prove (|38p we undertake the calculation of the Laplace transform of r. Observe first that ([55)) 
and ([36[) imply for 2p > d — 1 

m(2p,d,d)(dx) = (deta;)^'-"*^ [ ^ OFT^ ) IPdl^)^^:. (39) 
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In particular, in let us replace d by d — 1 and do 2p — d — 1. We get 



We can now write 



= -tr(st), 



idt) 



l,d- 


i,d- 


l){dx) . 


\ 


1 /n 




/ . 

/ 


e~ 


i (a/z 


j 








E 




^i) 


r{d/2) 


\k\1 

1 






E 




^i) 


r(<i/2) 


|«|! 

1 






E 




-i) 


r(d/2) 


\k\\ 

1 






E 




-i) 


r(d/2) 


\k\\ 

1 




7ri/2 


E 

Ke£d- 




-i) 


r(d/2) 


|«|! 

1 






E 

KSfd- 






T{d/2) 


\n\l 

1 




7ri/2 


E 




-i) 


r(d/2) 


|«|! 

1 




7ri/2 


E 


cMd- 


-i) 


r(d/2) 


1 





(40) 



tr (sux 



\K\\Td-i{n+^ 



-jry^-dx du (41) 



Vd-l{^+^) j 

{det[usu*]'^^)i <^ ^{[usu*Y['^)du 



«'«+|([usw*]i ^)du 



* n, (s"^)(dets-i)^ 



mi,...,md 



_„o(s-^)(dets-i)' 



, N-ii^ C'(k,o)(s ^) 



(42) 
(43) 
(44) 
(45) 
(46) 
(47) 
(48) 
(49) 



K,e£d- 



In this list equality (|4ip comes from , equality (I42p comes from (pSj) and equalities and 
(|T7| come from Lemma 4.3. In identities following (H^ we have replaced k by (mi, . . . ,md-i) for 
clarity. Formulas (|^^ and (H^ come from Lemma 4.2. and comes from Lemma 4.1. The proof 
of (|48p is more involved and is a consequence of formula (iii) of Theorem XIV 3.1 of FK where we 
replace (d, r, A,/i) respectively by l,d and 

d — 1 d — 3 d — 3 d — 1 
X = (mi H — ,m2H — ,...,md-i — , — ), 



A* = (- 



4 ' ' 4 
d-l d-1 d-3 



4 ' 4 

d- 3, 



, mi 



The fact that /i is a permutation of A and the reference above imply pS)) . The proof of comes 
from p|^y^.^ CK(/d-i) = CK.oild) implied by formula (|24p . An important point for this is simply 
e{K) = £(k,0). Finally gl]) proves dSS]) and Proposition 4.4 itself. □ 
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5 Convolution lemmas in the cone Vd 

We give the proof of Lemma 5.1 below though it is certainly a known fact. 

Lemma 5.1: In a Euclidean space E of dimension d let us fix a linear subspace F of dimension 
n. We choose randomly a linear subspace G of dimension k < d—n with the uniform distribution, 
that is the unique distribution on G such that G <^ uG for ah u e 0(ci). Then Pr(GnF ^ {0}) = 0. 

Proof: It is enough to prove the lemma for E = F = {0} x M''"" and k = d—n. Let Zi,. . . ,Zn 
be i.i.d. random variables in following the standard Gaussian distribution N{0,ld). Let G be 

the random linear subspace of E generated by the vectors Zi, . . . , Zn- Since for all u € 0{d) we 
have {uZi, . . . , uZn) ~ {Zi, . . . , Z^) clearly G ~ uG and G has the uniform distribution. Introduce 
the matrix 

M = [Zi,. . . , Zn] = {Zij)l<i<d, l<j<n 

whose columns are the vectors Zi, . . . , Z„. Then xiZi + - ■ •+a;„Z„ = MX where X = (xi, . . . , a;„)*. 
Now GOF 7^ {0} if and only if there exists a non-zero X such that the n first elements of MX are 
zero. In other terms, considering the square matrix Mi of order n defined by Mi = {Zij)i<ij<n, 
we have that G C] F {0} if and only if there exists a non-zero X such that MiX = 0, which 
happens if and only if det Mi = 0. Since the ti^ entries of the matrix Mi are independent N{0, 1) 
variables, the event det Mi = has probability zero and this proves the lemma. □ 

In the sequel we denote by Sb the set oi x Vd with b = rankx — 0, . . . ,d. 

Lemma 5.2: Let F be a random variable in Sb and assume that uYu* ^ Y for all u in the 
orthogonal group 0{d). Let xq € <S'o. Then xq + Y is concentrated on Sa+b if a + 6 < d or on 
Sd = Vd a + b > d. Furthermore if xq € Sc, and ii xq + Y is concentrated on Sa+b with a + b < d 
then c = a. 

Remark: U a + b = d and if Xq + Y is concentrated on Sa+b = Sd, Xq could be on any Sc with 
a < c< k. 

Proof: Apply Lemma 5.1 to F = xq'R'^ and to G = YW^. Then dim(F + G) = a + b if a + b < d. 
This implies that rank (xq+Y) < a +6. To see that rank (xq+Y) = a+b suppose that {xo + Y)R''' ^ 
E = F + G. Let x'q and Y' be the restriction of the endomorphisms xq and Y to the linear space 
E. Since xq and Y are symmetric, this implies that x'qE = F and Y'E = G. Since {x'q + Y')E ^ E 
there exists v G E \ {0} which is orthogonal to (a;Q + Y')E and this implies that {x'q + Y')v = 0. 
Since x'qV G F and Y'v G G and since F n G = {0} this implies that x'qV = Y'v = 0, and v is in 
Kcr(.TQ) n Ker(y ). But since F(BG = E {a, direct sum, not necessarily an orthogonal one) wc have 
also Ker(xo) © Ker(y) = E which implies Ker(xo) n Ker(F') = {0}. Thus v = 0, a. contradiction. 
Finally {xq + Y)R'^ = E = F + G and rank {xo + Y) = a + b. 

If a + 6 > d then F contains a subspace of dimension d, — b and dim(F + G) = d. Suppose now 
that Xq G Sc, and that xq+Y is concentrated on Sa+b with a-|-6<A;. Ifc-|-6<(i then, by the first 
part of the lemma, a;o + 1^ is concentrated on Sc+b- But Sa+b = Sc+b implies c = a. l{c+b>k 
then, by the first part of the lemma again, xo + Y is concentrated on Sd- This is impossible since 

Sa+b 7^ Sd- □ 

Lemma 5.3: Let fi and v be positive measures on Vd such that v is concentrated on Sb and v is 
invariant by the transformations x uxu* , u G 0{d). Let a = 0, . . . , d — 2 such that a + b < d. 
Then u is concentrated on Sa+b if and only if fi is concentrated on Sa- Furthermore v is 
concentrated on Sd = Vd if H is concentrated on Sd-b- 

Proof: For G Sb consider the distribution Ky^^ {dy) on Sb of the random variable UyoU* where 
U is uniformly distributed on the orthogonal group 0{k). Let Db the set of diagonal elements yo of 
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Sh of the form yo = diag(Ai, . . . , A;,, 0, . . . , 0) such that Ai > A2 > . . . > > 0. Then there exists 
a unique positive measure uq on such that the following desintegration holds 

v(dy) = / vo(dyo)Ky„{dy). 

It follows that 

{li*v){dx) = / iy{dy)fi{dx - y) = / n{dx - y) / V(i{dy(i)Ky^{dy) 

L 



/ Mdyo) / l^{dx - y)Ky^{dy) 



Therefore the measure ^ * Ky„ is concentrated on Sa+b for 1^0 almost all yo S -Dfc. From Lemma 
5.2 this implies that // is concentrated on Sa- 

<^= If /i is concentrated on Sa with a + 6 < A: it is an easy consequence of Lemma 5.2 that * 1/ 
is concentrated on Sa+b- O 

Lemma 5.4: Let a,b £ {!,..., d— 1} such that a + b < d. If m{a,a + b, d) exists, it is concentrated 
on Sa- 

Proof: From the Laplace transforms we know that m{a, a + b,d) * m{b, 0, d) = m{a + b,a + b,d). 
From Proposition 2.2 we know that m(o + 6, a + 6, rf) is concentrated on Sa+b- Since the Laplace 

transform of m{b,0,d) is (dets)"^/^ we know that m(b,0,d) is invariant by the transformations 
X i-> uxu* for any u G 0{d). By Lemma 5.3 we deduce that m{a, a + b,d) is concentrated on if 
it exists. □ 

6 m{d — 2,d — l,d) and m{d — 2, d, d) do not exist for d>2> 

In this section we prove Propositions 2.4 and 2.5. 

Proof of Proposition 2.4. Suppose that m{d — 2,d,d) exists. By Lemma 5.4 the measure 
m{d — 2, d, d) is concentrated on Sd,-2- By Lemma 5.3 the convolution 

•m{d — 2, d, d) * m(l, 0, d) = m(d — 1, d, d) 

is concentrated on Sd-i- This contradicts Proposition 4.4. □ 

Proof of Proposition 2.5. Suppose that m{d — 2,d — l,d) exists. By Lemma 5.4 the measure 
m{d — 2,d — l,d) is concentrated on Sd-2- Therefore there exists a positive measure m{dy) = 
mldyi, dyd-2) on M<*('^-2) = R<* x . . . x R'^ such that for all seVdwe have 



/ Q-{y'isvi+-+v2_^sva-2)^(^fj_y^ 

J^d(d-2) 



1 



,tr(s-i/(d-l,d)) 



(50) 



/Md(d-2) dets(''-2)/2 

We write more conveniently the elements y = {yi, . . . , yd-2) with the help of the transposed matrix 

'd 

= [Cl, . . • ,Cd] . 

• • • yd-2,d 

With this notation introduce the Gram matrix 

G{c) = G{ci,. ..,Cd) = {{cj, Ck))i<j,k<d 



y* = {yi,j) with d-2 rows yt,... 


yd-2 S'lid d columns Ci, . . . , 


Cd 




y*i 




2/1,1 ••• yi,d 




* 

y = 












. yd-2 . 




yd-2,1 ■ ■ ■ yd-2,d 





13 



and denote by m{dc) what we denoted by m{dy) before. We get 



e-tr(.G(c))„^(^^) = l^e*'^(^"^('^-i-'^)). (51) 

dets^" •^'/^ 

Equality (|5ip means that m(d — 2, d — 1, d){dx) is the image of m{dc) hy c^-^ x — G{c). 

Now in (|5ip we choose s — diag(l, si) where si is a symmetric positive definite matrix of order 
d — 1. We also desintegrate m{dc) by introducing a probability kernel i4r(c2, . . . , c^; dci) and a 
positive measure mi{dc2, . . . , dcd) such that 

e""''^" m{dci,dc2, ■ ■ ■ ,dcd) = mi{dc2, ■ ■ ■ ,dcd)K{c2, ■ ...Cd] dci) 
With these notations we can write 



e 

;d-i)(d-2) 



;d-i)(d-2) 



since if is a probability kernel. The last equality says that the image of mi {dc2 , ■ ■ ■ , dcd) by the map 
(c2, . . . , Cd) ^ X — G(c2, . . . , Cd) is nothing but m{d— 2, 1, d— l){dx). Denote G2 — G{c2, . . . , Cd) 
for simplicity. Since C2,...,Cd are vectors of a Euclidean space of dimension d—2 the rank 
of G2 is less than or equal to d — 2. To prove this elementary fact of linear algebra we use 
G2 G Vd-i- This implies that if a; = {x2, ■ ■ ■ ,Xd)* then 6*22; = if and only if x*G2X — 0. Since 
x*G2X = II X)i=2 ^i^ilp the linear space of x G M.'^^^ such that X)i=2 ^i'^i ~ ^ least dimension 

1, the kernel of the endomorphism of R'^^^ with matrix G2 has at least dimension 1 and its image 
has at most dimension d—2. This contradicts Proposition 4.4 which says that m{d— 2, d— 1, d— 1) 
has an absolutely continuous part and therefore charges matrices with rank d ~ 1. □ 
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